Affine LPV systems: realization theory, input-output equations and 
relationship with linear switched systems 
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Abstract — We formulate a Kalman-style realization theory 
for discrete-time affine LPV systems. By an affine LPV system 
we mean an LPV system whose matrices are affine functions 
of the scheduling parameter. In this paper we characterize 
those input-output behaviors which exactly correspond to affine 
LPV systems. In addition, we characterize minimal affine LPV 
systems which realize a given input-output behavior. Further- 
more, we explain the relationship between Markov-parameters, 
Hankel-matrices, existence of an affine LPV realization and 
minimality. The results are derived by reducing the problem 
to the realization problem for linear switched systems. In this 
way, as a secondary contribution, we formally demonstrate the 
close relationship between LPV systems and linear switched 
systems. In addition we show that an input-output map has a 
realization by an affine LPV system if and only if it satisfies 
certain types of input-output equations. 

I. Introduction 

The paper presents a Kalman-style realization theory for 
discrete-time affine LPV systems. An affine LPV system 
(abbreviated by ALPV) is linear parameter-varying systems 
whose matrices are affine functions of the scheduling pa- 
rameters. By the input-output behavior of an ALPV we will 
mean the input-output map induced by the zero initial state. 
The paper aims at answering the following questions. 

• How can we characterize those input-output maps which 
can be described ALPVs ? What is the role of Hankel- 
matrices in this characterization ? 

> What can be said about minimal ALPVs realizing the 
given input-output map ? What is the relationship be- 
tween minimal ALPVs, and reachability and observabil- 
ity of such systems ? Are all minimal ALPV realizations 
of the same input-output map isomorphic ? 

• How can we characterize the input-output equations 
solutions of which correspond to input-output maps of 
ALPVs ? 

In this paper we will show the following. 

• We prove that reachability and observability of ALPVs 
is equivalent to minimality and that minimal realizations 
of the same input-output map are isomorphic. Note that 
isomorphism in this setting means a linear state-space 
transformation which does not depend on the scheduling 
parameter. 



We also show that any ALPV can be transformed into 
a minimal one while preserving its input-output map. 
In addition, we characterize reachability and observabil- 
ity in terms of rank conditions for extended reachability 
and observability matrices. 

• We define the Markov-parameters as functions of the 
input-output map. We then show that the Hankel-matrix 
constructed from the Markov-parameters has a finite 
rank if and only if the corresponding input-output 
map has a realization by an ALPV. We show that the 
Kalman-Ho algorithm of [31] can be used to compute 
an ALPV realization from the Hankel-matrix, and we 
provide a bound on the size of the Hankel sub-matrices 
which guarantees correctness of the algorithm. 

• We also present a class of input-output equations which 
characterize ALPVs precisely: an input-output map is a 
solution of such an input-output equation if and only if 
it admits a realization by an ALPV . 

• Finally, as a secondary result, we establish a for- 
mal equivalence between the realization problems for 
ALPVs and for linear switched systems. The solution 
of the latter problem is known [23], [22], [21] and it is 
equivalent to that of recognizable formal power series 
and state-affine systems [7], [28], [12]. We then use 
realization theory of linear switched systems to derive 
a Kalman-style realization theory for ALPVs. 

Note that in this paper we consider ALPVs with a fixed 
initial state. Just as in the linear switched case [21], [22] it 
is possible to extend these results to the case of an arbitrary 
set of initial states. 

Motivation and novelty To the best of our knowledge, the 
paper is new. Many of the concepts (Hankel-matrix, Markov- 
parameters, extended reachability/observability matrix, etc.) 
used in the paper have already appeared before. However, 
what is truly novel in this paper is that it formulates a 
Kalman-style realization theory for ALPVs, while using 
the existing concepts from the literature. In addition, the 
equivalence between ALPV realizations and input-output 
equations is also new, to the best of our knowledge. 

A Kalman-like realization theory offers several benefits 
for system identification. It allows the characterization of 
identifiability and equivalence of state-space representations. 



The latter is important for model validation. Kalman-like 
realization theory also provides a tool for finding identifiable 
canonical parameterizations and characterizing the manifold 
structure of systems, including hybrid and nonlinear systems, 
[26], [29], [13], [20], [14], [15], [19]. In turn, this knowledge 
could be used for deriving new parametric identification 
algorithms, see [20], [14] for the linear case. Realization 
theory also leads to model reduction techniques, such as 
balanced truncation and moment matching [1]. This is also 
true for linear switched systems [25], [24] and ALPVs [31]. 

Finally, the paper formulates the precise relationship 
between the realization problems for ALPVs and linear 
switched systems. While this relationship is part of the 
folklore, it has not been stated formally yet. 

Relationship with existing work The field of identifica- 
tion of LPV systems is a mature one with a vast literature 
and several applications, without claiming completeness, we 
mention [36], [38], [35], [37], [18], [34], [30], [32], [16], 
[9], [8], [6], [3], [11]. As it was mentioned before, many 
of the concepts used in this paper were published before. 
In particular, the idea of Hankel-matrix appeared in [31], 
[36], [38], [35], [37]. However, [31], [36], [38], [35], [37] 
focuses on the identification problem, which is related to, 
but different from the realization problem studied in this 
paper. The Markov-parameters were already described in 
[31], [35]. In contrast to the existing work, in this paper the 
Markov-parameters and Hankel-matrix are defined directly 
for input-output maps, without assuming the existence of a 
finite dimensional ALPV realization. In fact, the finite rank 
of the Hankel-matrix represents the necessary and sufficient 
condition for the existence of an ALPV realization. The 
Kalman-Ho realization algorithm was discussed in [31], 
but it was formulated with the assumption that an ALPV 
realization exists. Moreover, the conditions under which the 
algorithm yields a true realization of the input-output map 
were not discussed in detail in [31]. Extended observability 
and reachability matrices were presented in [35], [31]. How- 
ever, their system-theoretic interpretation and relationship 
with minimality were not explored. 

Realization theory of more general linear parameter- 
varying systems was already developed in [30]. In [30] the 
system matrices are allowed to depend on the scheduling pa- 
rameter in a non-linear way. Moreover, in [30] no conditions 
involving the rank of the Hankel-matrix were formulated 
for the existence of a state-space realization. Hence, the 
results of [30] do not always imply the ones presented in 
this paper. The minimality conditions of [30] imply those 
of this paper. However, an ALPV may be minimal in the 
sense of this paper, and may fail to be minimal in the sense 
of [30]. Intuitively this is not at all surprising, since it is 
conceivable that by allowing more complicated dependence 
on the scheduling parameter we can get rid of some states. 

In particular, minimal ALPVs in the sense of this paper 
are related by constant state-space isomorphism. This is in 
contrast to [30], where the isomorphism relating state-space 
representations may depend on the scheduling parameter. 
Note that a minimal ALPVs in the sense of this paper need 



not be minimal in the sense of [30]. Hence, there might exist 
several state-space isomorphisms between ALPVs which 
are minimal in the sense of this paper. Some of these 
isomorphisms might depend on the scheduling parameters. 
However, the results of this paper imply that there will be 
a constant state-space isomorphism. This is also consistent 
with [17]. 

Although realization theory of ALPVs is quite similar 
to that of linear switched systems, there are important 
differences. In particular, there exist no parallel for linear 
switched systems of the equivalence between realizability 
and existence of input-output equations. In fact, ALPVs seem 
to behave more like state-affine systems [27], [28] for which 
an analogous result exists. 

It is well known that there is a correspondence between 
LPVs and LFT representations [36], [33]. In [2], [5], [4] 
the theory of recognizable formal power series was used to 
develop realization theory for LFT representations. In this 
paper we reduce the realization problem of ALPVs to that 
of for linear switched systems. The latter problem can also be 
solved by using recognizable formal power series [21], [22], 
[23]. Hence, there is an analogy between our approach and 
that of [2], [5], [4]. Note that the transformations between 
ALPVs and LFT representations involve non-trivial transfor- 
mations of the system matrices. Moreover, the resulting class 
of LFT representations seem to differ from the one in [2], 
[5], [4]. For this reason, it is unclear how the results of this 
paper could be derived directly from [2], [5], [4] and whether 
such an approach would be simpler than the current one. 

Outline In SflU we review the definition of ALPVs and 
the related system-theoretic concepts. In ^TTTl we establish 
the formal relationship between ALPVs and linear switched 
systems. In 31 VI we present a Kalman-style realization theory 
for ALPVs. Finally, in §|V] we present the input-output 
equations describing the behavior of ALPVs. 

Notation Denote by N the set of natural numbers including 
0. The notation described below is standard in automata 
theory, see [10]. Consider a (possibly infinite) set X. Denote 
by X + the set of finite non-empty sequences of elements 
of X, i.e. each w G X + is of the form w = a\ai ■ • ■ a&, 
ai,<X2, . . . ,afe G X, k > 0. The length of the sequence w 
above is denoted by \w\. We denote by wv the concatenation 
of the sequences w,v G X + , i.e. if w = a\ ■ ■ ■ aj. and 
v = v\---vi, a±, . . . , afe, vi, . . . , v\ G X, then wv = 
a\ ■ ■ ■ a^v i We denote by e the empty sequence. We 

define X* = X + U{e} as the set of all finite sequences of el- 
ements of X, including the empty sequence. By convention, 
|e| = 0, and the concatenation is extended to X* as follows: 
for all w G X*, we = ew = w. For each j = 1, . . . ,m, 
ej is the jth unit vector of R m , i.e. ej = (Sij, . . . ,S n ,j)> 
5j,j is the Kronecker symbol. If Z is a subset of a vector 
space, then SpanZ denotes the vector space spanned by the 
elements of Z. 

II. Discrete-time LPV systems 

In this section we present the formal definition of ALPVs 
along with a number of relevant system-theoretic concepts 



for ALPVs. 

Definition 1: A discrete-time affine linear parameter- 
varying system (abbreviated by ALPV) is of the form 

f x(t + 1) = EL(VW+VWKW m 
1 y(t) = Z%(c q x(t)) P M 

Here V C R D is the space of scheduling parameters, D 
is a positive integer, p(t) = (pi(t), . . . ,po(t)) & V is the 
scheduling signal, u(t) G R m is the input, y(t) G R r is 
the output and A q G R" x ™, B q G R" xm , C 9 G R rx ™, 
g G Q = {!> • • • j -D} are the system matrices. We will use 
the following short notation. 

(r > m > n,P,{(A 9> B <r ,C,)}^ 1 ) 
Notation 1: In the sequel, Q = {1, . . . , D}. 
The definition above also allows for affine dependence on 
the scheduling parameters. To this end, choose V to be of 
the form V = {(pi, ■ • ■ ,Pr>) | Pi = 1, (P2, ■ • ■ ,£>£>) 6 ^} 
for some set V C R- 0-1 . Moreover, if the affine hull of P 
equals R 15-1 , then the linear span of V will be equal to M. D . 
The latter property is important, because in the sequel we 
often use the technical assumption that V contains a basis 
of R D . 

Note that in our definition the output y t at time t does 
not depend on the input at time t. This restriction is made 
in order to simplify notation and most of the results can be 
easily extended to include direct dependence of y t on u t . 

Throughout the section, E denotes an ALPV of the form 
(QJ- The dynamics of E is driven by the inputs {u(t)}^ 
and the scheduling parameters {p(t)}^L . The state of the 
system at time t is x(t). If V = {e\, . . . , eu}, where 
denotes the ith standard basis vector, i = 1, . . . , D, then the 
ALPV E can be viewed as a linear switched system with the 
set of discrete modes being equal to Q = {1, . . . , D}. 

In order to enable formal discussion, we define a number 
of standard concepts such as input-output maps, reachability, 
etc. for ALPVs. 

Notation 2 (Generalized inputs): Denote U = V x R m . 
We denote by U* (resp. U + ) the set of all finite (resp. non- 
empty and finite) sequences of elements of U. A sequence 

w = 0(0), u(0)) • • ■ (p(t), u(t)) G U + , t > (2) 

describes the scenario, when the scheduling parameter p(i) 
and the input u(i) are fed to E at time i, for i = 0, . . . , t. 

Definition 2 (State and output): Let x G M. n be a state 
of S, Define the input-to-state map x^. x '■ U + — > M. n and 
input-output map yy,, x : U + — > W of E as follows. For any 
u> G W + of the form (O, define xs.x(w) as the state x(t) of 
E at time i, and define ys :X (w) as the output y(i) of E at 
time t, if the initial state a;(0) of E equals x, and the inputs 
{ u {i)Yi=o an d tne scheduling signal {p(i)}* =0 are fed to E. 
Note that for t = 0, .ts.£c(w) = x. 

The definition above implies that the potential input-output 
behavior of an ALPV can be formalized as a map 

/ :U + -> R r . (3) 



The value /(id) for u> of the form <j2j represents the output 
of the underlying black-box system at time t, if the inputs 
{u(i)Y i=0 and the scheduling parameters {p(i)}l =0 are fed 
to the system. This black-box system may or may not admit 
a description by a ALPV. Next, we define when an ALPV 
describes (realizes) /. 

Definition 3 (Realization): The ALPV E of the form (HJ 
is a realization of an input-output map / of the form (0, 
if / equals the input-output map of E which corresponds 
to the zero initial state, i.e. / = ys.o- The map ys,o will 
be referred to as the input-output map of E and it will be 
denoted by y-£. 

Similarly to [22], [21], the results of this paper could be 
extended to families of input-output maps and multiple initial 
states. However, in order to keep the notation simple, we deal 
only with the case when the initial state is zero. 

Definition 4 (Input-output equivalence): Two ALPVs Ei 
and E2 are said to be input-output equivalent, if y^ 1 = y-% 2 . 

Definition 5 (Reachability): Let E be an ALPV of the 
form ([TJi. We say that E is reachable, if the linear span of 
all the states of E which are reachable from the zero initial 
state yields the whole space R™. 

Definition 6 (Observability): The ALPV E is called ob- 
servable if for any two states Xi,X2 G R™, y^, Xl = Uy,. X2 
implies x\ = X2- 

That is, observability means that if we pick any two distinct 
states of the system, then for some input and scheduling 
signal, the resulting outputs will be different. 

Note that the concepts of reachability and observability 
presented above are strongly related to extended controlla- 
bility and observability matrices from subspace identification 
of ALPVs [35]. Later on, we will show that the ALPV is 
reachable if and only if the extended controllability matrix 
is full rank, and the ALPV is observable if and only if the 
extended observability matrix is full rank. 

Finally, we recall the notion of isomorphism for ALPVs. 

Definition 7 (ALPV isomorphism): Consider a ALPV Ei 
of the form (U) and a ALPV E 2 of the form 

Z 2 = (r,m,n a ,P,{(A a q ,B a q ,C a q )}° =1 ) 

with n a = n. A nonsingular matrix S G R" x ™ is said to be 
an ALPV isomorphism from Ei to E2, if 

Vq=l,...,D:A%S = SA q , B a q = SB q , C a q S = C q . 
Note that in the definition of an ALPV isomorphism, the 
state-space transformation S does not depend on the schedul- 
ing parameter. Finally, below we define what we mean by 
the dimension minimality of a ALPV. 

Definition 8 (Dimension): The dimension of E, denoted 
by dim E, is the dimension n of its state-space. 

Definition 9 (Minimality): Let / be an input-output map. 
An ALPV E is a minimal realization of /, if E is a 
realization of /, and for any ALPV E which is a realization 
of /, dim E < dim E. We say that E is minimal, if E is a 
minimal realization of its own input-output map y^. 



III. Relationship between linear switched 

SYSTEMS AND ALPVS 

In this section we establish a formal relationship between 
ALPVs and linear switched systems. We start by stating the 
following assumption. 

Assumption 1: In the rest of the paper, unless stated 
otherwise, we will assume that the linear span of elements 
of V equals R D , i.e. V does not belong to any of the proper 
linear subspaces of MP . 

Note that the assumption above is not restrictive. Indeed, if 
V belongs to a D dimensional proper linear subspace X of 
M D , then we can define a linear map S : MP — > U. D such 
that S is injective on X and replace the set of scheduling 
parameters by V = S(V). Since S is linear, the parameters of 
the resulting new LPV system will depend on the parameters 
in an affine way. 

Next, we introduce the concept of generalized convolution 
representation for input-output maps. This concept will allow 
us to concentrate on input-output maps for which there is a 
hope that they can be realized by ALPVs. 

Notation 3: Let p = p(0)---p(t) be a sequence of 
scheduling parameters and let v = qo ■ ■ ■ q t G Q + , qo ■ ■ ■ qt £ 
Q. Thenf = Pqo (0) Pqi (l)--- Pqt (t). 

Definition 10 (Convolution representation): Let / be an 
input-output map of the form (0. The map / has a gen- 
eralized convolution representation (abbreviated as GCR), 
if there exists a map S f : {v G Q + \v\ > 1} ->■ R rxm 
such that for each w G U + of the form (O, 
t-l 

/(«>)=£{ E Sf(v)f k:t }u(k), (4) 

k=0 VEQ+ ,|«|=t-fc+l 

where p k t = p(k)p(k + 1) ■ • -p(t). 

The convolution representation states that f(w) is linear in 
control input and that it is a homogeneous polynomial of 
degree one in the scheduling parameters. The values of the 
map play the role of the coefficients of this polynomial. 
Note that the concept of GCR above is a special case of 
impulse response representation (IRR) in [30]. Note that 
since in the ALPVs of interest the output at time t does not 
depend on the input at time t, the summation in (01i goes only 
up to t — 1. Below we show that is uniquely determined 
by / and that the existence of a GCR implies that without 
loss of generality we can assume that V = M. D . 

Lemma 1: If / has a GCR, then the map is uniquely 
determined by /. Moreover, there exists a unique extension 
f ext of / to U+t, where U ext = (K D x M m ), such that f ext 
also admits a GCR and S f = S f '^. 

Proof: [Proof of Lemma [T| The fact that f ext exists 
relies on the fact that (|4|l is defined for any values of 
p(0) , . . . , p(t) G M. D , and by noticing that the right-hand 
side of (HJi is a sum of terms, each of which multilin- 
ear in p(i), ■ ■ . , p(t), i = 0,...,t. Recall that function 
g(z\, . . . , Zk) is multi-linear, if for each i = 1, . . . , k, if we 
fix zi, . . . , Zi-i, Zi+x, ...,Zk and we vary only Zi, then g is 
a linear function of z,. Then set f ext as the value of the 
right-hand side of (|4j. If the value of f(w) is known for 



p(0), ■ ■ ■ ,p(t) where p(0), . . . ,p(t) run through a bases of 
M. D , then these values uniquely determine the value of the 
right-hand side of ©, and thus f ext exists and it is unique 
and S^ cmt = S*. Finally, the uniqueness of »V follow by 
noticing that S f (v)u = f ex t((e qo , u)(e qi , 0) • ■ • (e qt , 0)) for 
v = qo ■ ■ ■ qt, Qo, ■ ■ ■ , qt G Q, u € K m . ■ 
In the sequel, we will restrict attention to input-output maps 
which admit a GCR. This is not a strong restriction, since 
the input-output maps of ALPVs always admits a GCR. 

Lemma 2: The ALPV E of the form (|T|i is a realization 
of an input-output map / if and only if / has a GCR and 
for all v = qo ■ • ■ qt £ Q + , qo, ■ ■ ■ , qt £ Q, t > 

S (v) = C qt A qt _ 1 A qt _ 2 ■ ■ -A qi B qg . (5) 

If t = 1, then A qt l A qt 2 ■ ■ ■ A qi is interpreted as the 
identity matrix. 

Now we are ready to state the relationship between ALPVs 
and linear switched systems. To this end, we introduce the 
following notation. 

Notation 4 (Switched generalized inputs): Denote V sw = 
{d, ...,e D } and U sw = (V sw x R m ). 
Recall that we can view linear switched systems as a subclass 
of ALPVs, such that the space of scheduling parameters 
equals V sw . Potential input-output maps of linear switched 
systems are maps of the form / : U^ w — > R r such 
that / admits a GCR. Linear switched systems and their 
input-output maps in the sense of [23] correspond to linear 
switched systems and their input-output maps in the above 
sense, if one identifies the scheduling parameter e q with 
the discrete mode q G Q. We refer the reader to [23] for 
the notion of realization, minimality, observability, span- 
reachability, isomorphism. Alternatively, all these notions 
are special cases of the corresponding concepts for ALPVs, 
if one identifies linear switched systems as a subclass of 
ALPVs . Note that the concept of span-reachability from 
[23] corresponds to the concept of reachability as defined in 
Definition 

Definition 11: For each / : U + — > W admitting a GCR, 
define the associated switched input-output map 3(f) ■ 
U^w as follows. Let f^t be the extension of / to U^ xt 

as described in Lemma [Hand define 3(f) as the restriction 

of f ext to U+ w C U+ Xt . 

By noticing that = 5^ cxt = S ^- 1 we can in fact conclude 
that the correspondence between / and 3(f) is one-to-one. 
Next we will establish a correspondence between ALPVs 
and linear switched systems. 

Definition 12: Let E be a ALPV of the form (THi. 
Define the linear switched systems 6(E) associated 
with E as the linear switched system 6(E) = 
(r, m, n, V sw , {(A q , B q , C q )}^ =1 ). 

The following theorem collects the properties of the corre- 
spondence between linear switched systems and ALPVs. 
Theorem 1: 1) An ALPV E is a realization of the 

input-output map /, if and only if 6(E) is a realization 

of 3(f). 

2) For any ALPV E, dim 6(E) = dimE. 



3) Two ALPVs Si and £2 are isomorphic if and only if 
©(£1) is isomorphic to ©(£2). 

4) The ALPV £ is reachable, observable, minimal if and 
only if ©(£) is respectively reachable, observable, or 
minimal. 

Proof: [Sketch of the proof of Theorem [T) The only 
non-trivial statement is that © preserves reachability and 
observability. Let £ be an ALPV of the form (Q3. First we 
show that £ is reachable if and only if ©(£) is reachable. 
To this end, consider the map input-to-state map rr^.o : 
U + — > R" of £. Notice that sj^ can be extended to act 
on U^ xt and that for any input w £ U^ xt of the form 
©, xs.q is a sum of terms, each of which is multilinear 
in p(0), . . . ,p(t). Hence, the linear span of the values of 
x~£,o(w), w £ Utxt equals the linear span of values of 
££,o(w), w £ (Z x R m )+, where Z is a basis of R D . Since 
by Assumption Q] V contains such a basis of R-° and V sw 
is a basis of R. D , it follows that the linear span of xy,,o(w), 
w £ U + equals the linear span of Xe,o(w), w £ U+ w . Finally, 
notice that xs,q(w) — a;e(E),o(' u; ) f° r all w € Mtw Hence, 
£ is reachable if and only if ©(£) is reachable. 

Next, we show that £ is observable if and only if ©(£) is 
observable. To this end, notice that y-^ tX can be extended to 
U^ xt and that for any w £ U^ xt of the form ©, ys, x (w) is a 
sum of terms, each of which is multilinear in p(0), ■ ■ ■ ,p(t). 
Hence, ys yXl and yj2. X2 agree on U + , if they agree on any 
set (Z x R m )+, where Z is a basis of MP. Since V sw is 
a basis of R D and by Assumption [T] V contains a basis of 
R-°, it then follows that ys yXl and ys >X2 are equal on U + if 
and only if they are equal on Uf w . Notice that for i = 1,2, 
Ue(Z), Xi coincides with the restriction of ys. Xi to the set 
U+ w . This then implies that £ is observable if and only if 
©(£) is observable. ■ 

IV. KALM AN- STYLE REALIZATION THEORY 

In this section we exploit Section [TTT] and realization theory 
of linear switched systems [21], [22], [24], [23] to formulate 
a Kalman-style realization theory for ALPVs. 

We start with presenting a characterization of minimality. 

Theorem 2 (Minimality): An ALPV is minimal, if and 
only if it is reachable and observable. If two minimal ALPVs 
are equivalent, then they are isomorphic. 
The theorem above is a direct consequence of Theorem Q] 
and [23, Theorem 3]. 

Similarly to linear switched systems [23], one can con- 
struct example of an ALPV £ which is minimal (reach- 
able, observable), while none of the linear subsystems 
(A q , B q , C q ), q £ Q is minimal (resp. reachable, observable). 

Next, we present rank conditions for observability and 
reachability. To this end, recall from [31], [35] the defini- 
tion of extended reachability and observability matrices for 
ALPVs. That is, let £ be of the form (\}. We define the 
extended reachability matrices TZ{, i £ N for £ as follows: 
TZ Q = [B x , B 2 , . . . , B D ] and for all i £ N, let 

7^+i = [A x n h A, 11, 1,,/vV 



Similarly, we define the extended observability 
matrices Oi for £ recursively as follows: Co = 
[Cf , C 2 T , . . . , C T D \ T and for all i £ N, 



<+i 



[Aia 



T 

i ' 



Notice that lZ n -i equals the reachability matrix of the 
switched system ©(£) and 0„_i equals the observability 
matrix of ©(£). For the definition of reachability and 
observability matrices for linear switched systems see [23]. 
Hence, Theorem [TJ and [23, Theorem 4] yield the following 
rank conditions. 

Theorem 3: The ALPV £ is reachable if and only if 
rank7?. n _i = n, and £ is observable if and only if 
rank0 n _i = n. 

Theorem [3] yields algorithms for reachability, observability 
and minimality reduction of ALPVs. These algorithms are 
the same as those for linear switched systems [23]. 

Next, we present the necessary and sufficient conditions 
for the existence of a ALPV realization for an input-output 
map. To this end, we need the notion of the Hankel- 
matrix and Markov-parameters of an input-output map. In 
the sequel, / denotes a map of the form Q, and we assume 
that f has a GCR. 

Definition 13 (Markov-parameters): The Markov- 
parameter M*(v) of / indexed by the sequence v £ Q* is 
the following rD x Dm matrix 



M f (v) 



S f (lvl), 



Sf{Dvl) 
Sf(Dv2) 



(6) 



_Sf(lvD), Sf{DvD)_ 
That is, (v) can be viewed as a D x D block matrix, such 
that the (£,j)th entry of M*(v) equals S^(jvi), j,i £ Q. 

If / has an ALPV realization £, then from Lemma [2] it 
follows that Mf(v) can be expressed as product of matrices 
of £: if £ is as in ©, then M / (e) = CB and for all v = 
qi, . . . ,q k £ Q, k > 0, 



Mf(v) = CA qk A qk _ 1 .--A qi B, 



(7) 



where C=[Cf, C^]\ B = [B ly ...,B D ]. 

Note that the values of the map S* , and hence the Markov- 
parameters {M f(v)} v£ Q* can be obtained from the values of 
/. A naive way to compute 5^ is to compute the derivatives 
of / with respect to the scheduling parameter. It is easy to 
see that the Markov-parameters / and 3(f) coincide, i.e. 
Mf(v) = M d( -f\v), v £ Q*. Moreover, when applied 
to linear switched systems, the Markov-parameters from 
Definition [13] coincide with the ones in [23, Definition 12]. 

Note that the definition of Markov-parameters does not 
assume the existence of an ALPV realization of /. In fact, 
even if / does not admit a finite dimensional state realization, 
its Markov-parameters remain well-defined. The reason for 
this choice is that we want to use the Markov-parameters 
to characterize the existence of a finite dimensional ALPV 
realization of /. This will be achieved by constructing a 
Hankel-matrix from the Markov-parameters and by proving 



that / has an ALPV realization if and only if the rank of that 
Hankel-matrix is finite. Of course, for this to make sense, we 
have to define the Markov-parameters and the Hankel-matrix 
as objects which are well-defined even in the absence of a 
finite dimensional state-space representation. 

In order to define the Hankel-matrix of /, we will intro- 
duce a lexicographic ordering on the set Q*. 

Definition 14 (Lexicographic ordering): Recall that Q = 
{1, . . . , D}. We define a lexicographic ordering -< on Q* as 
follows. For any v, s S Q*, v -< s holds if either (a) \v\ < \s\, 
or (b) < \v\ = \s\ = k, v s and the following holds: 
v = qx - ■ - q k , s = si - ■ ■ s k , q x , . . . ,q k , s±, . . . , s k G Q, and 
for some / £ {1, . . . , k}, qi < s; with the usual ordering of 
integers and qi = Sj for i = 1, . . . , I — 1. Note that -< is a 
complete ordering and 



= {vi,v 2 , ■ • ■} 



(8) 



with vi -< i>2 -< • • ■■ Note that v\ — e and for all i e N, 

Definition 15 (Hankel-matrix): Define the Hankel-matrix 
Hf of / as the following infinite matrix 



Hi 



Mf{vwi), M f (v 2 v{), M f (v k vi), 
Mf(yiV2), Mf(v 2 v 2 ), Mf(v k v 2 ), 
Mf( Vl v 3 ) Mf(v 2 v 3 ), Mf(v k v 3 ), 



i.e. the rD x mD block of Hf in the block row i and block 
column j equals the Markov-parameter M*(s), where the 
word s = VjVi € Q* is the concatenation of the words Vj 
and Vi from ©. 

Note that H s = H 3 ^ and the definition of the Hankel- 
matrix coincides with the one for linear switched systems 
[23, Definition 13]. 

Theorem 4 (Main result on existence): The map / has a 
realization by an ALPV if and only if / has a GCR and 
rankff/ < +oo. Any minimal ALPV realization of / has 
dimension equal to rankify. 

The theorem above is a direct consequence of Theorem Q] 
and [23, Theorem 5]. 

Finally, we prove the correctness of the Kalman-Ho- 
like realization algorithm for ALPVs from [31]. A similar 
algorithm was formulated for linear switched systems in 
[24], [23]. To this end, we need the following definition. 
For every LeN, denote by N(L) = J2j= VJ the number 
such all the sequences v £ Q* of length at most L. Due 
to the properties of lexicographic ordering, it follows that 
{«!,..., %(L)} = {v G Q* I \v\ < L}. 

Definition 16: Denote by Hf j^^j the N(L)rD x 
N(M)mD upper-left sub-matrix of Hf. 
If / is realized by an ALPV E, then H ftLM = G L K M , 
where Ol is the Kth extended observability matrix and TZm 
is the Mth extended reachability matrix of E. In this case 
H/.L.M coincides with the Hankel-matrix defined in [31]. 
The Kalman-Ho algorithm goes as follows. Compute the 
factorization 

Hf : L,L+i = OR 



such O € W DN ^ xn , R £ K»xm£>N(L+i) and ran kO = 
rankR = n for n = recnkH fx.L+i- One way to compute 
this factorization is by SVD decomposition as in [31], i.e. if 
Hfx.L+i = USV T is the SVD decomposition of Hf^,L+i 
where S is the diagonal part, then set O = US 1 ! 2 and R = 
S 1/2 V T . Let R be the matrix formed by the first N(L)mD 
columns of R. For each q 6 Q, let R g be the n x N(L)mD 
matrix, such that the jth n x mD block column of R g equals 
to the fcth n x mD block column of R, where k is such that 
vjq = v k . Here v k and Vj are the jth and kth elements of 
the lexicographic ordering ([SJ. Construct E of the form ([]]) 
such that [Bi, B]j\ equals the first mD columns 



ofR [Cj 



^ 2 1 



of O and A q = R 9 R + 



Cj] equals the first rD rows 



where R is the Moore-Penrose 
pseudoinverse of R. 

Theorem 5: If rankiJ/ ^ i = rankif/, then E computed 
by the algorithm above is a minimal realization of /. The 
condition rankiJ/ ^ ^ = rankify holds, if there exists an 
ALPV realization E of / such that dimE < L + 1. 
The theorem above is a direct consequence of Theorem Q] 
and [23, Theorem 6]. 

V. Input-output equations for ALPVs 

In this section we use the results of realization theory 
to establish a relationship between ALPVs and input-output 
equations. In the sequel, / is assumed to be an input-output 
map f : U + — > R r and it is assumed that f admits a GCR. In 
order to avoid excessive notation, in this section we assume 
that r = 1. However, all the results can easily be extended 
to several outputs. 

Definition 17 (Input-output equations): An ajfine polyno- 
mial equation E(P, Y, U) of order n is a polynomial in 
variables P = ]/'.,}, „. r Q . Y = {rj™ =0 , U = 



j /<=!,. ..,n,3=l, 



such that 



E(P, Y, U) = J2 Qi(P)Yj + E E LuPWu (9) 



3=0 



i=l j=l 



where Qq(P), Qj(P), Lij(P) are polynomials, i = 
1, . . , , n, j = 1, . . . ,m and Qo(P) ^ 0. 

Definition 18: Assume that E is an affine polynomial 
equation of the form (|9). Then the input-output map / is 
said to satisfy the equation E, if for each w of the form (O 
wit t > n, E(f, w) = 0, where E(f, w) denotes the value of 
i7(P, Y, U) with the following substitution Pij = pj(t — i), 

U hl = u,(t-i) Y l = /((p(0),u(0)) •••(?(*- *),«(*-*)) 
for j G Q, I = 1, . . . , m, i = 0, . . . , n. 

Theorem 6: Assume that the set of scheduling parameters 
V is an open subset of R D . The input-output map / has a 
realization by an ALPV if and only if / satisfies an affine 
polynomial equation of the form (O. 

In [30] it was shown that input-output maps of LPV systems 
with a meromorphic dependence on parameters correspond 
to input-output maps which satisfy linear autoregressive 
equations with respect to outputs and inputs. The coeffi- 
cients of these autoregressive equations were meromorphic 
functions of the time-shifted scheduling parameters. Affine 



polynomial input-output equations represent a special case 
of the autoregressive equations of [30]. Theorem [6] says that 
input-output maps described by these type of equations (and 
which, in addition, admit a GCR) correspond precisely to 
input-output maps realizable by ALPVs. 

The proof of Theorem [6] is an adaptation of the proof of 
the analogous statement for state-affine systems [27], [28]. 
The proof is divided into several lemmas, proofs of which 
are presented in the appendix. 

Lemma 3: If the interior of V not empty, then / satisfies 
the input-output equation (0 if and only if its extension f ext 
from Lemma Q] satisfies (0. 

From Lemma |3] it follows that without loss of generality, we 
can assume V = M. D . 

Assumption 2: In the sequel, we assume that V = M. D . 

For any sequence p = p±p 2 ■ ■ -Pk & T y+ , pi, . . . ,pk £ 
k > define the map /£ : U + — > K as follows: 

V»6W+: fP{w) =/(w(pi,0)(p2,0)---(pfc,0)) 

Recall that w(pi,0) ■ ■ ■ (pk,0) denotes the concatenation of 
the sequence w with the sequence (pi, 0) • • ■ (pk, 0). Intu- 
itively, f-(w) equals the response of /, if first we feed in 
the inputs and scheduling parameters prescribed by w and 
then for the last k time steps we feed in the zero input and 
the scheduling parameters pi,...,Pk- 

Lemma 4: There exists an affine polynomial input-output 
equation E of the form (0 such that / satisfies E, if and 
only if there exists polynomials Qj(P), i = 0, . . . ,n such 
that Q 7^ 0, and for any p\, . . . ,p„+i e V, 

n 

J^Q j (pi,---,Pn+l)f PlP2 "' Pn+ ^ (10) 
3=0 

Before formulating the next statement, recall the set of all 
maps g : U + —> W forms a vector space with respect to 
point-wise addition and multiplication by scalar. 

Lemma 5: The map / satisfies ([TUt for some Qj, j = 
0, . . . ,n if and only if W f = Span{/£ | p £ (R D )+} is 
finite dimensional. 

Lemma 6: The input-output map / has a realization by a 
ALPV if and only if Wf = Span{/£ | p e (R D ) + } is finite 
dimensional. 

The proof of Lemma [6] boils down to showing that there 
is a linear isomorphism between Wf and the linear space 
spanned by the rows of the Hankel-matrix Hf of /. Hence, 
Wf is finite dimensional if and only if rankff/ < +oo. By 
Theorem HI the latter is equivalent to the existence of an 
ALPV realization of /. Theorem [6] follows from the lemmas 
above as follows. From Lemma [6j / has a realization by a 
ALPV if and only if Wf is finite dimensional. By Lemma [5] 
and Lemma|4] the latter is equivalent to existence of an affine 
polynomial equation of the form (0 such that / satisfies E. 

VI. Conclusion 

We have presented realization theory for the class of affine 
LPV systems. In addition, we have shown that realization 
theory of this class of LPV systems is equivalent to that 
of for linear switched systems. We have also presented 



an equivalent input-output representation for affine LPV 
systems. 

Acknowledgements We thank Roland Toth for the nu- 
merous remarks and suggestions. 
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Appendix 

Proof: [Proof of Lemma [3) Note that E(f ex t,w) is 
a polynomial expression in p(0), u(0)), . . . , (p(t),u(t)) £ 
U ex t for each w of the form (f2]). Since E(f ex t,w) = for 
all w € U + and U is an open subset of U ex t, it then follows 
that E(f ext , w) = for all to € U ext . ■ 

Proof: [Proof of Lemma |4) By substitution it is clear 
that if (0 holds, then ( TTOb holds. Conversely assume that 
( fTOl l holds. With the notation of ©, notice that for F* = 
/((p(0), u(0)) • • • (p(t - i), «(* - »))• * = 0, . . . , n, 

y. — j?p(t-n)p(t-n+l) — y p(t-i) (<y)_|_ 
n 

+ R j,Mt-j),...,p(t-i))u(t-j), 

3=i+l 

where v = (p(0), tt(0)) • • • (p(t — n— 1), u(t — n— 1)), and 
Rj,i, j = n, . . . , i + 1 are suitable polynomials. Consider 
the expression J2™=o <9i(p(* -n),... ,p(t))Yi. By grouping 
together the terms (Q%Rj i)(jp{t — n), . . . ,p(t))u(t — j) in a 
suitable way, we can obtain polynomials Ljj,, j = 1, . . . , n 
and I = 1, . . . , m such that then (0 holds. ■ 
Proof: [Proof of Lemma[5) Assume that / satisfies ( TTOb . 
Since Qo ^ 0, there exists an open an dense subset ZCP 
such that for any pi , . . . ,p n +i £ Z, Q (pi, ■ ■ -,p n +l) i= 0. 
Since Z is open, there exists bi,...,brj £ Z such that 
bi,...,b D span V = WL D . Let B = {b u . . . , b D }. Then for 



any pi,... ,p n +i £ B, Qo(pi, • ■ -,Pn+x) 0- By dividing 
CGD by Q (p 1 ,...,p n+1 ), it follows that is a 

linear combination of fPi'"Pi > i = 1, ...,n. It then follows 
that for any p\, . . . ,pk+ n +i £ B, fPv"P*+n+y is a linear 
combination of f Pl , i = 1, ...,n. Hence, we then 
obtain that f-, p £ B + belongs to the linear span of the 
finite set {ft \ ~p £ B+, \p\ < n}. But /pi-** is multi- 
linear in pi, . . . ,pk, since / has a GCR. Hence, W/ is the 
linear span of the finite set {ft | p g B + , \p\ < n}. Hence, 
W/ is finite dimensional. 

Conversely, assume that Wf is finite dimensional. For 
each v = q 1 ■ ■ ■ q k £ Q + , q 1 ,...,q k £ Q, denote by 
/" the map / e «i "' > e ik . As it was noted above, /P 1 *"'?* is 
multilinear in pi,...,pk and hence W/ equals the linear 
span of f Zl , . . . , f Zd for some z\, . . . , Zd £ Q + . Notice that 
for any p £ (R D )+, \p\ = fc, ft = Z ve Q+M=k /V- 
Since for every v £ Q + , f" is a linear combination of 
f Zi , i = l,...,d, there exist polynomials P^fc is kD 
variables, such that /pi-p* = p j,k(j>i,- ■ -,Pk)f Z] for 

any p 1: . . . ,p k £ M. D . 

Consider now the rf x (d + 1) polynomial matrix D^+i 
in variables Xi = (X^i, . . . , X^d), i — 1, 2, . . . , d + 1 such 
that (i, j)the entry of D^+i equals Pij(Xi, . . . ,Xj), Let's 
view Dd+i as a matrix with elements in R(Xi, . . . , Xd+i). 
Here, M.(Xi, . . . , Xd+\) is the quotient field of the polyno- 
mial ring M.[Xi, . . . , Xd+i}- Since D^+i has only d rows 
and d + 1 columns, the columns of D^+i must be lin- 
early dependent. It then follows that there exist polynomials 
Dj,Nj £ R[X 1 ,...,X d+ i}, Nj jt 0, j = l,...,k*, such 
that D k * 7^ and Pid~N~ — 0- multiplying the 

equation above by the product of N\ ■ ■ ■ Nk* we get that 

k* 

yi = l,...,d:J2 P i,j R j=° (ID 

j'=i 

for some polynomial i?i, . . . i?^*, i?^* 7^ 0. Notice that the 
polynomial Pij depend only on the variables X\, . . . ,Xj, 
hence i?i , . . . , Rk* can be chosen to be polynomials 
only in Xx,...,X k *. If k* = 1, then P ifl = and 
hence f p = for all p £ V. Hence, fP^-^( w ) = 
f Pk (w(p 1 ,0)---(p k - 1 ,0)) = for for all to £ W+ 
Pi,...,Pk, k > 0. Then ( fTOl l holds for n = 1 with any 
choice of Qi and Qo- If fc* > 1, then set n = k*, 
Q l = Rk'-i, i = 1, • • • k* - 1. Using the fact that /pi-p< = 
Ej=i ■ • ■ >ft)/ Zj and (HB, it then follows that ([TOj 

holds for all p u ... ,p k £ R D . ■ 
Proof: [Proof of Lemma |6) Denote by H the linear 
span of the rows of the Hankel-matrix Hf . Notice that each 
element of TL can be viewed as a sequence of 1 x Dm 
matrices. We define the linear map $ : Wf — > H as follows: 
$(/?_) = (if^ ; iJ t , 2 ,...), such that for each v £ Q* , 

H v = Y, If' ,---,P sD ]Mf(vs). 
seQ+,|s|=|p|-i 

In other words, H v = [H Vt i , H V; d], where H VlQ = 



SseQ+ |s|=|p| S^(qvs)p s . Moreover, for any w of form (f2]), 

t-i 

Hence, it is clear that $ is an injective linear map. Moreover, 
the row of Hf indexed by the integer I = (i — l)D+q,q G Q, 
1 = 1,... equals $(/ e «) if i = 1, or <&(f e n '" e "k < e «), if i > 1 
and gi, . . . , Qfe e Q are such that Vi = qi ■ ■ ■ qu, where Vi is 
ith sequence of the lexicographic ordering ([8}. Hence, $ is 
a linear isomorphism from W/ to the space spanned by the 
rows of Hf. The rest of follows from Theorem [4] 



